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ABSTRACT 

The central star-forming region in a blue compact dwarf galaxy, II Zw 40, was observed in 
the 340 GHz (880 urn) band at ~ 5 arcsec (250 pc) resolution with the Submillimetre Array 
(SMA). A source associated with the central star-forming complex was detected with a flux of 
13.6 ± 2.0 mJy. The structure is more extended than the beam in the east-west direction. The 
SMA 880 [am flux is analyzed by using theoretical models of radio spectral energy distribution 
along with centimetre interferometric measurements in the literature. We find (i) that the SMA 
880 (J.m flux is dominated (~ 75 per cent) by free-free emission from the central compact star- 
forming region, and (ii) that the contribution from dust emission to the SMA 880 [J.m flux is at 
most 4 ± 2.5 mJy. We also utilize our models to derive the radio-FIR relation of the II Zw 40 
centre, suggesting that free-free absorption at low frequencies (v < several GHz; A > several 
cm) and spatial extent of dust affect the radio-FIR relation. 

Key words: dust, extinction — galaxies: dwarf — galaxies: evolution — galaxies: individual 
(II Zw 40) — H II regions — submillimetre: galaxies 



1 INTRODUCTION 

Blue compact dwarf galaxies (BCDs) generally host compact and 
ongoing star formation activities in metal-poor and gas-rich en- 
vironments dSargent & Searlel 1 19701 : Ivan Zee. Skillman. & Salzen 
1998). Because both low metallicity and rich gas content indicate 
an early evolutionary stage, BCDs can be used as nearby labo- 
ratories of primeval galaxies which should exist at high redshift. 
In some BCDs, the most active class of star formation is tak- 
ing place in super star clusters (SSCs) (Turner, Ho, & Beckl l 19981 : 
IKobulnickv & Johnsonll 19991: 1 Johnson & Kobulnickv 2003), which 
are expected to 'mimic' the starburst in high-redshift primeval 
galaxies. 

II Zw 40 is a well studi ed BCD with a low o xygen abun- 
dance, 12+log(0/H) = 8.13 JThuan & Izotovl2005l) . This galaxy 
hosts a high star formation activity associated with the SSCs in 
the centre. At centimetre wavelengths, the central star-forming re- 
gion in II Zw 40 is compact and optically thick for free-free ab- 
sorption , and is categorized as a 'supemebula' or 'ultradense H I I 
region' JTurner. Ho. & Becld[i998l : IKobulnickv & Johnsonll 19991) . 
Such a dense and compact star for mation activity is called 'ac- 
tive' mode in iHunt et al.1 J2003I) and Iffirashita & Huntl d2004l) . II 
Zw 40 is also classified as a Wolf-Rayet galaxy: the Wolf-Rayet 
feature in dicates that the typic a l age of the current starburst is a 
few M yr dVacca & ContHl 1992b . fBuckalew. Kobulnickv. & DufouJ 
(2005) derived an age of 2.6 Myr from the Ha and H/3 equiva- 
lent widths. The Br7 equivalent width also shows age < 3 Myr 
dVanzi et al.l2008f) . Stellar spectral synthetic models support young 
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ages ~ 2 Myr, although there is an underlying old stellar popula- 
tion dWestera et aT]|2004l) . 

In order to trace dense 'embedded' star-formation activ- 
ities, optically thin star formation indicators are usefu l, such 
as far-infrared (FIR) dust lumi nosity (e.g. iKennicuttl 1 19981 : 
llnoue. Hirashita. & Ka mava 2000) and radio (thermal plus non- 
thermal) luminosity dCondonlll9920 . Indeed, there is a correlation 
between FIR and radio luminosities in nearby star-forming galax- 
ies (e.g., de Jong et al. 1985; Helou, Soifer, & Rowan-Robinson 
1985), which is naturally explained if bo th luminos ities are strongly 
connected with star formation activities ( Volk 1989). However, the 
studies o f radio-FIR re lation are biased to objects with FIR detec- 
tion (e.g. Condon 1992), which means that a significant dust enrich- 
ment has already occurred. Therefore, the evolution of radio-FIR 
relation in young primeval galaxies is not yet clear. The gas den- 
sity, the magnetic field strength, and the energy density of cosmic 
ray ele ctrons affect the evolution of radio emission on various time- 
scales dHelou & B icav 19931: iNiklas & Beckj[l997l : lMurphvll2009l: 
lLacki & Thompsonll20100 . whi le the dust enrichment p lays a role 
to increase the FIR luminosity dHirashita & Hunt 2008). 

There have been some observational studies on the evo- 
lution of radio-FIR relation along the cosmic age. The radio- 
FIR relation at moderate and high redshifts (z < 4) is 
broad l y similar to that at t he local Univer s e (e.g. JG arrett 
20021: iGruppioni et alj 120031: llbar et alj 120081 : iMurphyj 1200*91 



_ _ _ JL. _ 

Micha lowski. Watson. & Hi orth 2010 ), although there is also a 

slight indication of evolution 1 Vlahakis, Eales. & Dunnd 120071 : 



ISevmour et al. 2009; Michalowski et al.ll20101) . The disadvantage 
of high-z observations is the difficulty in deriving quantities re- 
lated to galaxy evolution (age, metallicity, etc.) with high ac- 
curacy. On the other hand, studies of nearby metal-poor dwarf 
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galaxies as 'laboratories' of primeval galaxies provide an al- 
ternative way to approach the properties of pr imeval g alaxies. 
iHopkins. Schulte-Ladbeck. & Drozdovskvl J2002h and IWu et al.l 
(2008) conclude that the radio-FIR relation of star- forming dwarf 
galaxies is similar to that of normal galaxies in spite of the dif- 
ference injnetallicity (see a l so|Klein. Weiland. & Brinks! 19911) . al- 
though |Cannon_et_alJ (2005, 2006) find a significant deviation from 
the canonical radio-FIR relation for some individual dwarf galax- 
ies. 

In this paper, we do not take a statistical way, but investi- 
gate a single object in details. II Zw 40 is suitable for studying 
radio-FIR emission in a metal-poor young object, since as men- 
tioned above it hosts an extremely young (< 3 Myr) active star 
formation with a metallicity of 1/ 4 Zq (for the solar metallicity, 
we adopt 12 + log O/H) = 8.69; lLoddersll2003h . In order to spot 
the young star-forming component, a high resolution is necessary. 
The star-forming region is resolved well at centimetre wavelengths, 
while there is no information on FIR-submillimetre (submm) dust 
emission on such a small scale. A submm interferometric observa- 
tion is desired for the purpose of resolving the dust emission in the 
star-forming region. Therefore, we observed II Zw 40 by the Sub- 
millimetre Array (SMA; Ho, Moran, & Lo 2004), and we report on 
this observation in this paper. This observation, combined with cen- 
timetre radio interferometric data in the literature, will enable us to 
obtain the radio-FIR emission properties (or radio-FIR relation) in 
the very young star-forming region with a low metallicity. 

This paper is organized as follows. We explain the observa- 
tions and the data reduction in Section[2] and describe observational 
results in Section[3] In Section[4] we interpret the results along with 
the radio data in the literature by using theoretical models. In Sec- 
tion[5] we analyze the contribution from dust emission to the SMA 
flux. We discuss the radio-FIR relation based on the SMA observa- 
tion in Section[6] Finally, we conclude in SectionfT] The distance to 
II Zw 40 is assumed to be D = 10.5 Mpc (cz — 789 km s _1 with 
Ho = 75 km s _1 Mpc -1 ). At this distance, 1 arcsec corresponds 
to 50.9 pc. 



2 OBSERVATIONS AND DATA 

The SMA observation of II Zw 40 was carried out in the 340 GHz 
(880 pin) band on 2010 March 30 in the subcompact configuration. 
Six antennas were used with projected antenna separations between 
9.5 and 25 m. The receivers have two sidebands, the lower and up- 
per sidebands, which covered the frequency ranges from 331.0 to 
335.0 GHz, and from 342.9 to 346.9 GHz, respectively. The visibil- 
ity data were calibrated with the MIR package, with Mars as a flux 
calibrator (with an adopted flux of 1194 Jy), quasars J0423— 013 
and J0730— 116 as amplitude and phase calibrators, and quasar 
3C273 as a band pass calibrator. The calibrated visibility data were 
imaged and CLEANed with the MIRIAD package. The synthesized 
beam has a size of FWHM 5.2" x 4.4" (265 pc x 224 pc) with a 
major axis position angle of ~ 60°. The beam roughly matches the 
size of the re gion of active star fo rmation in the Hubble Space Tele- 
scope i mag e (Calzetti et al. 2007) and the Very Large Array (VLA) 
image l lUlvestad, Johnson, & Nefnl2007l) . The largest angular scale 
sampled by this observation is 19". 

In order to estimate the contribution of free-free emission to 
the 880 pm SMA flux, we use radio data in the literature (Ta- 
ble QJ. The single-dish observations measure the fluxes from the 
entire galaxy, while the interferometric data are sensitive to angular 
sizes smaller than # max . We also adopt the 'matched beam' VLA 



Table 1. Radio data used in this paper. 



V 


n a 

p raax 


Flux 


Image b 


Fit c 


Ref. d 


(GHz) 


(arcsec) 


(rajy) 








1.4 


— 


30.0 ±0.5 


S 




1 


1.5 


— 


29.9 ±0.7 


S 




2 


1.5 


— 


30.5 ± 1.5 


S 




3 


4.8 


— 


22 ±3 


S 




4 


5.0 


— 


22 ±4 


S 




1 


5.0 


10 


15 ± 1 


I 




5 


5.0 


4 


9± 1.5 


M 


V 


5 


5.0 


— 


21.0 ± 1.9 


S 




2 


8.3 


7 


12 ± 1 


I 




5 


8.3 


4 


10 ± 1.5 


M 


V 


5 


15 


4 


14 ± 1.5 


I 


V 


5 


25 


— 


18 ±4 


S 




4 



a The maximum size scale that is well sampled by the interferometric 
observation. 

'S', T, and 'M' represent single dish flux, interferometric flux, and 
matched beam flux, respectively. 
c If the data is used to fit the radio emission models in Section [4] this 

column is marked with V. 

d Re fere nces. 1) Ijaffe, Perola & Tarenghil Il978l); 2) iKleinetalJ 
il99lh: 3)lDeeg et alj Jl993l) : 4) lKlein. Wielebinski. & Thuanl Jl984l) : 
5) lBeck et all |2002|) . 



Table 2. FIR and submm data used in this paper. 



A 

(unri 



Flux 
(mJy) 



Image ' 



Ref. 1 



60 
65 
70 
90 
100 
140 
160 
160 
450 
450 
850 
850 



6610 ± 700 

6900 ± 700 

5580 ± 280 

6600 ± 700 

5800 ± 900 

3700 ± 700 

3400 ± 900 

3140 ± 430 

400 ± 90 

248 ± 81 

90 ±10 

98 ±14 

13.6 ±2.0 



a 'S' and T represent single dish flux and 
interferometric flux, re spectively. 
b References . 1) iHunt et alj J2005h; 2) 
Hiras hitaet all I200SI); 3) lEngelbracht et all 
<2008l) ; 4) lGalliano et all 120051) ; 5) this paper. 



fluxes obtained with (u, v) data restricted to baselines greater than 
20kA (i.e. sensitive to structures smaller than 4 arcsec) jBeck et all 
2002). These matched fluxes are used to constrain the properties 
of the central active star formation. The discrepancy between the 
matched data and the single-dish data is larger at lower frequen- 
cies, which implies that the spectral slope is steeper in the diffuse 
medium than in the central star-forming region. 

For the information on dust emission, we also use FIR data 
in the literature as listed in Table [2] where the fluxes except for 
our SMA data are all measured by single-dish observations and 
represent those of the entire galaxy. We only consider the 'large 
grain' component which is in radiative equilibrium with the ambi- 
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Figure 1. SMA 880 |j.m continuum image of II Zw 40. Solid contours are 
2, 3, 4, 5, and 6 <r (1 ct = 1.15 mJy beam -1 ), while dotted contours are 
— 2 and —3 a. The beam is shown in the lower right corner. 



ent interstellar radiation field (see Section [5Tt . Thus, we do not use 
the data at mid-infrared and shorter wavelength s, where very small 
grains and PAHs d ominate dust emission (e.g. iDraine & Lill200ll : 
iDopita et al . 2006) and the spectral energy distribution (SED) is 
sensitive to grain size distribution. 



3 RESULTS 

Figure Q] shows the obtained image. The source is clearly detected 
with a peak intensity of 7.4 a (1 <7 = 1.15 mJy beam - ) and a 
total flux of 13.6 ± 2.0 mJy. There is an extended feature in the 
east-west direction with FWHM = 7.1 arcsec (the width of the 2 a 
contour is ~ 10 arcsec). 

Submm emission in galaxies is usually dominated by dust 
thermal radiation and contaminated by free-free emission. In order 
to estimate the contribution from free-free emission, we use the 
VLA 2 cm continuum data, which is sensitive to spatial scales up 
to 4 arcsec (Beck et al. 2002). We adopt the emission on this scale 
as the contribution from the central star-forming region associated 
with the SSCs. At this wavelength, the single-dis h flux is 18.5 mJy , 
while the VLA flux is 14 ± 1.5 mJy (Table [Tl lleck et alj|2002n ; 
thus, 76 per cent of the total radiation comes from the central star- 
forming region at 2 cm. Since the emission at 2 cm (14 ± 1.5 mJy) 
is dominated by free-free emission (Beck et al. 2002), we extrap- 
olate the fl ux to 880 ^m by assuming a frequency dependence of 
oc u~ 01 dOsterb rock 1989) and estimate the contribution of free- 
free emission to the 880 \xm flux to be 10 ± 1 mJy. Thus, 74^} 
per cent of the flux detected at 880 |j.m by SMA is free-free emis- 
sion from the central star-forming region. The contribution from 
free-free emission is further investigated by theoretical models in 
Section [4] The residual (4 ± 2.5 mJy) is possibly dust emission 
or diffuse (more extended than 4" but less than ~ 10") free-free 
emission. 



4 THEORETICAL MODELS 

In order to characterize the radio continuum radiation from a star- 
forming region, thermal free-free radiation from H II regions, non- 
thermal synchrotron emission from supernova remnants (SNRs), 
and dust thermal emission should generally be taken into ac- 
count. Here, we adopt theoretical models for thermal free-free 
and no n-thermal synchrotron components from iHirashita & Hunt! 
(2006) to interpret the emission from the central star-forming re- 
gion in II Zw 40. Below we briefly summarize their models. We 
adopt the same values for the physical parameters as those in 
IHirashita & Hunt! ( 120061) unless otherwise stated. We also add a 
simple calculation for dust emission in Section|5] 



4.1 Basic setups for the star formation 

The star formation is assumed to occur in a uniform star-forming 
region, whose initial number density of hydrogen nuclei is uho- 
The mass finally converted into stars, Mo, is also given. The star 
formation rate (SFR) as a function of time, ip{t), is written as 



ip(t) = (esFM /Jii)i 



't/ta 



, where t is the time measured 



from the onset of the star formation (thus, we assume ip = for 
t < 0), and iff is the free-fall time-scale estimated b y tiho (i.e. 
the S FR is regulated by the dynamical time-scale; e.g. lElmegreenl 
2000). We simply introduce an efficiency factor esf, which means 
that the gas is converted into stars on a time-scale of ta/esF — 
1.4 x 10 6 (e S F/0.1)~ 1 (n HO /10 5 cm -3 ) -1/2 yr. We treat n Hu and 
Mo as independent parameters. Throughout this paper, we assume 
a Salpeter initial mass function (IMF) with a stellar mass range of 
0.1-100 M . 

Since star formation does not proceed on a single free-fall 
time-scale because of some regulatio n factors such as magnetic 
fields, turbulence, radiation, etc . (e.g. [P rice & Bate 2009), we ten- 
tatively adopt esF = 0.1 (e.g. lLada. Lombardi. & Alvesl l2010)rl 
Although esF may be uncertain, the quantities actually constrained 
by the observational data are the gas mass converted into stars and 
the gas density (see Section |4~3t . As long as these quantities are 
focused on, the results are insensitive to esf- Moreover, ta/esF is 
shorter than the lifetimes of massive stars for most of the parameter 
range in this paper; in this case, the star formation is regarded as 
instantaneous as long as we consider the evolution at i > tg/esF- 
If the star formation is instantaneous, total stellar mass Mo, not 
the detailed star formation history, determines the total luminos- 
ity. Special remarks will be made when we discuss ages too young 
(< 1 Myr) to regard the star formation as instantaneous. 

4.2 Radio SED 

Thermal free-free radiation is modeled by considering the evolu- 
tion of the size and density of H II region. In this paper, the re- 
combination coefficient excluding captures to the ground level i s 
assumed to be a (2) = 1.73 x 10" 13 cm 3 s" 1 (e.g. lspitzerlll978h . 
We assume the density in the H II region to be uniform, and neglect 
the effects of dust in the H II region. For details on the assumptions, 
we refer to Hirashita & Hunt (2006). 

We consider the time evolution of the hydrogen number den- 
sity in the Hll region, iih, and the radius of the ionized region, 

The efficiency factor egp is different from the star formation efficiency, 
which is usually defined as the fraction of the finally formed stellar mass to 
the total initial gas mass. In Lada et al. (2010)'s notation, egp corresponds 
tol//~ 0.18. 
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Figure 2. Evolutions of the hydrogen number density njr (Panel a), the ionized radius r\ (Panel b), the number of ionizing photons emitted per unit time JVi on 



(Panel c), and the cumulative number of SNe V (Panel d). Solid, dotted, and dashed lines are for rajjo 
M = 3 x 10 6 M . 



3 X 10 , 10 J , and 3 X lCr cm d , respectively, with 



n, according to the evolution of the number of ionizing photons 
emitted per unit time, JV; on . We consider the decrease of nu by 
pressure-driven expansion, which is assumed to continue until n 
starts to decrease by the death of massive stars. In calculating the 
SED of thermal free-free component, the gas te mperature in the 
HIT region is assumed to be T — 13,000 K ( Thuan & Izotov 
2005). We also take free-free absorption into account. In particu- 
lar, for an H II region in such a dense region as the II Zw 40 centre, 
free-free absorption is important. The free-free optical depth rg i s 
estimated by using the emission measure EM as iHunt et alj|2004h 
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where EM is estimated in terms of the electron number density n c 
and the radius of ionized region n as EM = 4n e n/3. We assume 
that the absorbing medium is intermixed with the emitting material; 
that is, the escape fraction is assumed to be [1 — exp(— Ta)]/rg. 

The non-thermal radio luminosity of SNRs is estimated as 
^nt(i / ) = lntT n t{v/5 GHz)~°' 5 7(£), where j(t) is the rate of 



core-collapse supernovae (SNejj estimated by assuming that stars 
more massive than 8 Mq become SNe, and l n tT nt is the ra- 
dio energy emitted by a SNR over its entire lifetime (fluence) 
at 5 GHz. The fluence p ossibly depends on the ambient den- 
sity JArbutina & Urosevicl 120051) . We adopt Z n tT n t(5 GHz) = 
7.6 X 10 22 W Hz -1 yr, since this value fits the radio SED of 
SBS 0335—052, which also hosts dense and compact H II regions 
JHirashita & Hunt! 120061) . Free-free absorption is also applied for 
all the non-thermal component by a screen geometry: i.e. the ob- 
served non-thermal luminosity is L nt — L l l t e~ T{t . This mini- 
mizes the contribution from the non-thermal component especially 
at low frequencies. If the non-thermal component is not absorbed 
efficiently, the existence of the non-thermal component is more 
severely excluded for the II Zw 40 centre (see Section [43}. Syn- 
chrotron self-absorpti on is not imp ortant for SNRs in the frequency 
range of interest (e.g. ICondonll 19921) . 



In this paper, we only consider core-collapse SNe originating from mas- 
sive stars. 
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We also define the cumulative number of SNe, T, as 

T(i)= / j(t')dt'. (2) 

Jo 

r(t) is used as an indicator of the dust production in SNe (Section 

E3. 



4.3 Theoretical results for the radio SED 

We present the results calculated by the framework described 
above. We concentrate on the parameter ranges relevant for the cen- 
tral star formation in II Zw 40. As shown later, who ~ 10 5 cm -3 
and M ~ 3x 10 6 Mq fit the radio SED if such a young age 
~ 3 Myr as indicated by the optical observations (see Introduc- 
tion) is adopted. The radio luminosity is almost proportional to Mo, 
so the constraint on Mo is rather severe. Thus, we mainly investi- 
gate the case of Mo = 3x 10 6 Mq. For the density, we examine 
n H0 = 3 x 10 4 , 10 5 , and 3 x 10 5 cm -3 . 

In Fig. [2] we show the key quantities (riH, n, Mon, and T) 
as functions of time. The ionized region radius n monotonically 
increases and the density nu decreases until t ~ 3 Myr because 
of the pressure-driven expansion and the increase of Ni on . After 
that, the expansion stops because iN/ion decreases. The initial SFR is 
V>(0) = esFMo/tg ~ 2.2(£sf/0.1)(Mo/3 x 10 6 M )(n HO /3 x 
10 5 cm' 3 ) 1 ' 2 Mq yr" 1 . The SFR measured by the Ha lumi- 
nosity is roughly ~ 1 Mq yr" 1 jvan Zee et al.ll 1 9981 ; I Vanzi et al.l 
2008). Since the SFR exponentially decays on a time-scale of 
t ff /esF ~ 1.4(esF/0.1) _1 (n H o/10 5 cm" 3 )" 1 ' 2 Myrinourmod- 
els, we obtain the SFR averaged for 3 Myr as ~ 0.91 Mq yr" 1 



(tiho = 10 5 cm" 3 is assumed), which is near to the value ob- 
tained from the Ha line. As shown later, the SFR assumed also 
reproduc es the interferometric radio continuum flu x. For the stel- 
lar mass, iBuckalew, Kobulnickv. & Dufoun ( 120051) derived 6.3 x 
10 6 Mffl from a st ellar spectral synthesis model (Starburst 99; 
iLeitherer et alJI 19991) with a Salpeter I MF of stellar mass range 1- 
120 Mq and a metallicity of 1/5 solar. I Vanzi et alj J2008T) derived 
1.7 x 10 6 Mq with the same spectral synthesis model but with a 
Kroupa IMF. Our stellar mass (Mo ~ 3 x 10 6 Mq) is bracketed 
by those two results, which means that our stellar mass is consistent 
with those in the literature within the uncertainty in the IMF. The 
evolution of V is used for the discussion of dust production by SNe 
in Section 153] T rapidly increases after t — 3 Myr, when the first 
SNe occur. 

In Figure [3^, we show the radio SEDs at t — 1, 3, and 5 Myr 
with riHo = 10 5 cm" 3 and M = 3 x 10 6 Mq. As the age be- 
comes older, the peak shifts to lower frequencies, because the free- 
free optical depth becomes smaller as the H II region expands. At 
t — 1 and 3 Myr, the emission is completely dominated by free- 
free emission, and at t = 5 Myr, the synchrotron component be- 
gins to contribute to the emission and the spectrum slope changes. 
For comparison, the VLA 'matched' data whose (u, v) coverage 
is restricted to baselines greater than 20kA (i.e. sensitive to struc- 
tures smaller than 4 arcsec) are adopted (three triangles in Fig. [3} 
as representative fluxes from the central star-forming region. 

The density strongly affects the frequency at which the flux 
peaks because free-free absorption is sensitive to the density. In 
Fig.[3p, we show the SEDs at t = 3 Myr for various initial densities 
(tiho) with Mo = 3 x 10 6 Mq. We observe that the peak position 
of the SED is indeed sensitive to the density. The rising spectrum of 
the matched data (triangles) is consistent with free-free absorption. 

It is possible to search for the best-fit values of Mo and jiho 
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Figure 3. Radio SEDs. (a) Solid, dotted, and dashed lines show the SEDs 
at 1, 3, and 5 Myr, respectively. The initial hydrogen number density and 
the total mass converted into stars are fixed with njjo = 10 5 cm -3 and 
Mq = 3x 10 6 Mq, respectively. The observational data (Tableff} are also 
presented. The filled square represents our SMA measurement at 880 u.m 
for the c entral star-forming region. Diamonds s how single-dish m easure- 
ments by Ijaffe et alj h978Tl . iKlein et alj J1984I) . iKlein et alj h99ll). while 
asteris ks and triangles indicate the VLA interferometry data bv lBeck et alj 
12002). The triangles at 6, 3.6, and 2 cm are the matched VLA fluxes ob- 
tained by restricting (u, v) data to baselines greater than 20kA (i.e. sensitive 
to the structures smaller than 4 arcsec) IBeck et al . 2002). The asterisks at 6 
and 3.6 cm are the VLA fluxes synthesized by all the (u, v) data, which are 
sensitive to structures up to 10 arcsec and 7 arcsec, respectively, (b) Same 
as Panel (a) but for the SEDs at 3 Myr with various initial densities. Solid, 
dotted, and dashed lines are for njjo = 3 X 10 4 , 10 5 , and 3 X 10 s cm -3 , 
respectively (n H = 1.8 X 10 3 , 3.2 X 10 3 , and 5.4 X 10 3 cm" 3 at 3 Myr, 
respectively), (c) Same as Panel (a) but for the best-fit model parameters 
(Mo and nno) listed in Table [3] The fitting is applied for the three trian- 
gles (the matched VLA fluxes). The solid, dotted, and dashed lines show 
the best-fit solutions for t = 1, 3, and 5 Myr, respectively. 



for each age. The matched VLA data are adopted (three triangles 
in Fig.O for the \ 2 fitting (one degree of freedom), since our mod- 
els are applicable to the central star-forming region. The best-fit 
solutions are shown in Table [5] and Fig. [5J:. In the case of t — 1 
Myr, a large Mo is required since only 29 per cent of Mo is con- 
verted into stars at t — 1 Myr. As seen in Fig. [3};, the spectral slope 
at u > 15 GHz for t = 5 Myr is different from that at t < 3 
Myr because of the contribution from the non-thermal component. 
However, the spectral slope at v > 15 GHz in the central part of 
II Zw 40 is not well constrained because of the lack of the high- 
resolution data at millimetre wavelengths. For all the three ages, a 
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Table 3. Best-fit solutions. 



Age 

(Myr) 


M 
(M ) 


"HO 

(cm" 3 ) 


n H 
(cm" 3 ) 


X 2 


1 
3 

5 


1.1 x 10 7a 
3.5 x 10 6 
4.1 x 10 6 


2.2 x 10 4 
1.2 x 10 5 
1.2 x 10 5 


4.7 x 10 3 
3.7 x 10 3 
3.9 x 10 3 


2.6 

2.0 
2.9 



a Only 29% of M (i.e. 3.2 X 10 6 Mq) is converted into stars 
at t = 1 Myr. 

stellar mass of ~ 3-4 X 10 6 Mgjj and a current number density 
in the Hll region of 4-5 x 10 3 cm -3 are obtained for the best-fit 
values. 

If the age is about 3 Myr as sugge sted from the opt ical and 
near-infrared observations (Introduction; I Vanzi et alj|2008f) . 75 per 
cent of the 880 (im flux obtained in our observation is explained 
by free-free emission according to the best-fit SED at 3 Myr (the 
flux at 880 |J.m is 10.2 mjy in the model, while the observed flux 
is 13.6 mjy). The difference (~ 3.4 mjy) is likely to be due to the 
dust and/or diffuse (i.e. not associated with the compact H II region) 
free-free emission. The dust emission in the central star-forming 
region is modeled and examined in the next section. 



5 CONTRIBUTION FROM DUST IN THE CENTRE 

5.1 Dust emission associated with the central star-forming 
region 

Here we model the dust emission associated with the central star- 
forming region in II Zw 40. We do not model the entire system, so 
that the following FIR luminosity and dust temperature are valid 
for the central part of II Zw 40, and are not representative for the 
global quantities. We consider the wavelength range where large 
grains, which achieve radiative equilibrium with the ambient stellar 
radiation field, govern the emission (A > 100 ixm), since the large 
grains dominate the total dust mass (e.g. lGalliano et alj|20051) . We 
do not treat very small grains contributing to mid-infrared emission. 
For simplicity, we adopt a shell-like geometry for the 
dust distribution. Thi s simplification, which is also assumed in 
iTakeuchi et ail J2003L 120051) and iGalliano et all J2005h , helps to 
decrease the number of physical parameters that affect the emer- 
gent dust emission SED in relatively minor ways. If the dust is dis- 
tributed in a thin shell at a distance i?dust from the centre, the dust 
optical depth, Tdust, for the radiation from stars is estimated as 

3M dust ^ (3) 



"7~dus 



l&rJ^e 



where Mdust is the total dust mass in the shell, a is the grain radius 
(uniform spherical grains are assumed), and s is the grain mate- 
rial density. Note that -Rdust will typically be larger than n since 
the dust associated with surrounding neutra l or molecular g as con- 
tributes significantly to the FIR luminosity dXu et al. 19921) . In this 



paper , we adopt a — 0.1 vim and s — 3 g cm _J (e.g. lDraine & Led 
1 19841) . We estimated the absorption cross section of a grain with 
the geometrical one because the m ajor heating sourc e is ultraviolet 
(UV) radiation from OB stars (e.g. lBuat &XvJll996l) (i.e. the grain 
radius is comparable to the wavelength). 



3 At t = 1 Myr, only 29 per cent of Mq (i.e. 3.2 X 10 6 Mq) is converted 
into stars, while at t > 3 Myr, most of Mq is converted into stars (TablefJJ. 



We assume that the UV luminosity is equal to the bolometric 
luminosity of the OB stars. The OB star luminosity is estimated as 
a function of time: 



£oE>(t) 



*0 



dm / dt' L(m) (f>(m)ip(t — t'), 



(4) 



where L(m) is the main sequence luminosity of a star with mass m 
a nd 4>{m) i s the IM F. For L{rri), we adopt the zero metallicity case 
in lSchaererl (2002J) to consider a low-metallicity evolutionary stage, 
and L(m) may be 2 times smaller if we assume solar metallicity. 
We assume that the OB stars are located at the centre. This assump- 
tion is valid when we consider the central SSCs as the source of UV 
radiation in II Zw 40. Since the absorbed UV light is reprocessed 
into FIR wavelengths, the dust FIR luminosity is estimated by 



£fir = (1 — e 



')£ob. 



(5) 



For simplicity, we adopt a single-temperature approximation for the 
dust emission, so that the monochromatic luminosity is expressed 
as 



^Fm(^) = 47TK„ Mdust B„ (Tdust) 



(6) 



where B v is the Planck function, k„ is the mass absorption coeffi- 
cient of the dust, and Tdust is the dust temperature. We a ssume that 
re„ = 0.7(^/340 GHz) 2 cm 2 g _1 jjames et al.|[2002l) . The dust 
temperature is determined so that it satisfies the total luminosity 
constraint: 



£fir = 



Ltir{v)&v. 



By combining Equations l[5}-(|7), we obtain 



1-e" 



£oB/M dus t = 1.09 x 10- 5 T d 6 ust . 



(7) 



(8) 



The dust emission SED, Lfir(z/), is determined as follows. 
Rdust is given. Then, a certain test value of Md us t is assumed. By 
using Equation l[3j, Tdust is obtained. Consequently Td us t is derived 
from Equation (8). Note that the age has to be assumed to evaluate 
£ob (Equation [4)- Finally, the dust SED Lfir(z/) is obtained by 
Equation {5J. If the predicted flux LFnt(i/)/(47r.D 2 ) overproduces 
(underproduces) the observed flux, we decrease (increase) Mdust 
and repeat the above procedure. 



5.2 Comparison with observational data 



As mentioned in Section 1431 the difference between the predicted 
flux and the measured one by SMA at 880 \un should be dust emis- 
sion or diffuse free-free emission. To simplify the discussion, we 
interpret that all the difference is dust emission, which means that 
the following values for Afd us t and Lfir should be taken as up- 
per limits. In the SMA image (Fig. Q}, the extension of the 880 
(j.m emission is ~ 10 arcsec. That is, Rdust ^S 250 pc. Any 
more extended dust component contributes to the difference be- 
tween the SMA flux and the SCUBA flux. Also, R dus t < 60 pc 
is rejected because the dust temperature is so high that the 60 
(im flux exceeds the IRAS flux. Thus, we investigate -Rdust = 
100 and 250 pc as representative cases. For the free-free compo- 
nent, we adopt the best-fit model for t — 3 Myr (Table [3), i.e. 
n m = 1.2 x 10 5 cm" 3 and M = 3.5 x 10 6 Mq (the dotted 
line in Figure [3};). For the dust emission, we adopt the dust SED 
model in Section |5~fl and we adjust Mdust so that it reproduces 
the SMA observational data according to the formulation in Sec- 
tion O we obtain (Mdust, Tdust) = (1.8 x 10 4 Mq, 45 K) 
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Figure 4. FIR-submm SEDs. Solid, dotted, and dashed lines represent, re- 
spectively, the total, the free-free component, and the dust component cal- 
culated by the models for the central region in II Zw 40. For the free-free 
component, the best-fit SED to the matched VLA data for t = 3 Myr (i.e. 
the dotted line in Fig.[5J;) is adopted. The dust mass is adjusted to reproduce 
the SMA data point at 880 ^m (filled square). Two cases for the radius of 
dust distribution are examined: i?dust = 100 and 250 pc for the upper and 
lower solid lines, respectively. As upper limits for the flux in the central re- 
gion, we also show observational data for the entire II Zw 40 system (open 
diamon ds), which are tak en from Vader et alj 41993ft for 60 and lOOjjm 
(7R45 ), lHuntetalJ fe005h for 450, and 8 50 yim (SCUBA) iGalliano et at] 
l2005» for 450 and 8 50 p (SCUBA). lEngelbracht et alj (20081) for 70 
and 160 \im (Spitzer), Iffirashita et alj 120081) for 65. 90, 140, and 160 um 
(AKAR1). The observational data are summarized in Table|2] 



and (2.7 x 10 4 M , 35 K), for i? dus t = 100 and 250 pc, 
respectively. We present the SEDs in Figure [4] As upper lim- 
its for the flu x from the centra l part, w e also show the dat a 
taken b y IRAS JVader et alJll993h . AKARI JHirashita et al.ll2008l) 
Spitzer JEngelbracht et alj20 08) and SCUBA JGalliano et alj2005l ; 
iHunt. Bianchi. & Maio lino 2005) for the fluxes from the entire sys- 
tem. Although mid-infrared Spitzer data are available with resolu- 
tions comparable to our SMA data, mid-infrared emission comes 
from a different dust component, very small grains, as mentioned 
in Section |5~T| Thus, we only concentrate on FIR and submm wave- 
lengths where large grains dominate the dust emission. 

Fig.|4]shows that a significant fraction of the 60-100 urn flux 
(71 per cent of the IRAS 60 vim flux and 54 per cent of the IRAS 
100 um flux) is explained by the emission from the central star- 
forming region for i?dust = 100 pc. If we adopt i?dust = 250 pc, 
the contribution from the central star-forming region is smaller (22 
per cent of the IRAS 60 urn flux and 3 1 per cent of the IRAS 100 urn 
flux) because of a lower dust temperature. In any case, an additional 
emission from a diffuse dust component is necessary to explain the 
FIR fluxes of the entire system, especially, at wavelengths longer 
than 100 urn. 



5.3 The origin of dust in the centre 

Now we consider the origin of dust in the II Zw 40 cen- 
tre. One of the natural explanations for the dust mass derived 
above (~ 2-3 x 10 4 Mq) is production in previous episodes 
of star formation (i.e. preexisting). First, before investigating 
the possibility of preexisting dust, we examine how much dust 



can be produced in the current star formation activity. Because 
the age is young, SNe are the only source of dust grains and 
evolved low-mass stars such as asymptoti c giant branch (AGB) 
stars are negligible ( Valiante et al. 2009; G all. Andersen. & Hiorthl 
I2OIII) . Theoretical studies suggest that ~ 0.1-1 M^ of du st 
grains condense in a SN lIKo zasa. Hasegawa, & Nomoto 198 



Todini & Ferrara 2001; Nozawa et al. 2003; Bianchi & Schneide: 



ider 



2007; Nozawa et al. 2007, 2010). Infrared and submm observations 



of SNRs also detec t ed ~ 0.01-1 Mq o f dust (e.g..[Rho et al. 2009; 
lDunneetal.ll2009l ; iGomez et al.ll2009l ; iMatsuura et alj|201lh . al- 
though quantitative significance of SNe t o the dust enrichme nt in 
the early galaxy evolution is still debated JNozawa et alj|2010l . and 
references therein). If we adopt 1 Mq as the 'maximum' dust yield 
per SN, the total dust mass 1-3 x 10 4 Mq is explained only if at 
least 10 4 SNe occur. (If we consider dust destruction in SNe, more 
SNe are necessary to explain the total dust amount.) According to 
Fig-IU r does not reach this number even at t = 10 7 yr. Thus, we 
conclude that the excess at 880 urn observed in the centre of II Zw 
40, if the emission comes from dust, should be contributed from 
preexisting dust. 

As for the possibility of preexisting dust, we may be observ- 
ing the preexisting dust itsel f, or dust grains after the growth by the 
accretion of metals. iHirashita & Kuol a201 11) derive the dust growth 
time-scale for silicate (a similar time-scale is obtained for carbona- 
ceous dust): 



2.1 x 10 7 yr 



(a 3 )/(« 2 ) 
0.1 urn 



1Z 



© 



"n 



. 10 3 cm- 



- l (Te-.y 1 ' a (±y 1 (9) 



where (a 3 ) and (a 2 ) are the averages of a 3 and a 2 (a is 
the grain radius) for grain size distribution, Z is the metallic- 
ity, nu is the hydrogen number density, T gas is the gas tem- 
perature, and S is the sticking efficiency of the relevant metal 
species onto the dust surface. We assume (a 3 ) / (a 2 ) — 0.1 um, 
Z = 1/4 Zq (see In troductio n), n H = 10 5 c m" 3 (Tabled, 
y gaa = 50 K JWilson. Walker . & Thornievl 1 1 997h and S = 0.3 
JLeitch-Devlin & Willi ams 1985; iGrassi et aU201 II) . Then, we ob- 
tain Tg row ~ 8.4 x 10 5 yr, which is comparable to the star forma- 
tion time-scale. Therefore, if the preexisting grains survive before 
the current star formation episode, they can grow by the accretion 
of gas-phase metals. 



6 RADIO-FIR RELATION 

The 880 um emission in the II Zw 40 centre is dominated by free- 
free emission. At this wavelength, the emission is usually domi- 
nated by dust on a galactic scale. Thus, we expect that the properties 
of FIR dust emission and radio free-free emission in the II Zw 40 
centre are different from global galaxy properties. In this section, 
we examine this expectation quantitatively by using the results in 
the above sections. Hereafter, we use the term "FIR luminosity" for 
the dust emission luminosity integrated over 8-1000 um. 

In Fig. [5] we show the relation between the monochromatic 
luminosity at 1.4 GHz (Li.4GHz) and the FIR luminosity Lfir for 
the central star-forming region in II Zw 40. The radio luminosity at 
1.4 GHz is adopted from the best-fit model f or t = 3 Myr (Figure 
|3p), while the FIR luminosity is evaluated for i?dust = 100 and 
250 pc (Section [5.2| Fig. |4]l. Since all the difference between the 
model free-free flux and the observed 880 urn flux is assumed to 
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come from dust, the FIR luminosities estimated from the models 
should be taken as upper limits. 

For comparison, we also plot the observational data of BCDs 
in Fig. [5] The sample is taken from Hopkins et al. (2002) for the 
FIR and 1.4 GHz luminosities of the entire system (i.e. global lumi- 
nosities). We only adopt the objects with dete ctions by the NRAO 
VLA Sky Survey (NVSS) at v = 1.4 GHz JCondon et alii 19981) 
and by the IRAS at A = 60 /^m (large open squares). We also 
adopt the radio flux me asured by the Faint Images of t he Radio Sky 
at Twenty cm (FIRST; iBecker. White. Helfand[l993) (small open 
squares). The FIRST flux is syste matically smaller than the NVSS 
flux. The probable reason which iHopkins et alj (2002) suggest is 
the difference in their sensitivity to extended emission: FIRST im- 
ages tend to miss extended emission whose angular size is larger 
than about 2 arcmin. The FIR luminosity of the BCD sample is 
estimated based on the IRAS 60 /im and 100 /im fluxes obtained 
from NASA/IPAC Extragalactic Database (NED). As an observa- 
tional estimate of the F IR luminosity Lfir, we adopt an empiri- 
cally derived formula bv lNagata et al.l (2002), who estimate the to- 
tal dus t luminosity a t A ^ 4 u.m by using the IRAS 60 and 100 (jm 
fluxes, INagata et alj ( 120021) take into account the contribution from 
wavelengths longer than the IRAS bands by integrating the modi- 
fied blackbody spectrum, but they do not consider the contribution 
from A < 40 [im. Therefore, the FIR fluxes may be underesti- 
mated by < 30 per cent. However, INagata et alj J2002h 's method 
ca n avoid a sig nificant underestimate in a simple IRAS luminosity 
bv lHelou et alj Q988), who only considered the contribution from 
the luminosity in the IRAS FIR bands[jlf a sample BCD is not de- 
tected at 100 u.m, we utilize the upper limit at 100 [im to estimate 
an upper limit of Lfir. The data with upper limits are shown by 
crosses in Figure \5\ (Li. 4 ghz is taken from the NVSS data). We 
also show the global lu minosities of II Zw 40, whose 1 .4 GHz an d 
FIR data are taken from ljaffe et"al] dl978l) and lVader et ail d 19931) . 
respectively. 

In order to quantitatively discuss the possible deviation from 
the standard radio-FIR relation for BCDs, we adopt the radio-to- 
FIR ratio as usually used (e.g. lCondon 1992). Here we define 51.4 
and qi5 as 

Lfir 



gi.4 



Qis 



log 



= log 



3.75 x 10 12 W 

Lfir 
3.75 x 10 12 W 



log 



— log 



/L1.4 GHz\ 

HwhFtJ 

L15 GHz 



WHz" 



(10) 



(11) 



The average and the standard deviation of 91.4 are calculated for the 
sample detected by both IRAS (60 and 100 um ) and NVSS, whil e 
those of gig are determined for the sample in iHunt et alj (2005). 
The averages of gi.4 and gis are 2.39 and 2.89, respectively, and the 
standard deviations are 0.19 and 0.26, respectively. We also show 
the lines with qx.4 = 2.39 ± 0.19 and gi 5 = 2.89 ± 0.26 to show 
how much the radio-FIR relation in the central part of II Zw 40 is 
deviated from the reference radio-FIR relation for the BCDs. 

From Fig. [5] we observe that the central region of II Zw 40 has 
a suppressed 1.4 GHz luminosity and is located below the radio- 
FIR relation of global BCDs, if the FIR luminosity of II Zw 40 has 
a value near to the upper limit. This is interpreted to be due to free- 
free absorption in our models (Fig.O. Our SMA observation is not 
sensitive to the possible diffuse synchrotron radiation that would 
not be affected by free-free absorption. Indeed, the relatively high 



4 For II Zw 40, the global FIR luminosity estimated by Nagata et al 
2002l) 's method is 1.2 X 10 9 Lq, while that estimated by Helou et al 
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Figure 5. Radio-FIR relation for the central part of II Zw 40 (large filled 
circles) in terms of the global relations for BCDs. Radio frequencies of 1.4 
and 15 GHz are adopted in Panels (a) and (b), respectively. Two cases for 
the radius of dust distribution are presented for Lfir: i?dust = 100 and 
250 pc for the right and left points (corresponding the two SEDs in Figure 
|4)- The arrow indicates that the FIR luminosities estimated from the models 
are taken as upper limits since all the excess over the model free-free emis- 
sion is assumed to be dust emission. The observational data for the global 
emission from BCDs in Panel (a) are taken from Honkin s et alj (2002J) 
(large/small open squares are for NVSS/FIRST data and cros ses are for 
non-detection at 100 p.m), while those in Panel (b) are taken from lHunt et alj 
120051) . The entire II Z w 40 system is also plotted (1.4 GH z, 1 5 GHz, and 
FIR d ata are taken from ljaffe et alll978tlKlein et al.ll99l and lVader etaT] 
1 19931 respectively). We also show the lines with 91.4 = 2.39 ± 0.19 in 
Panel (a) and 515 = 2.89 ± 0.26 in Panel (b) (solid lines with dotted lines 
for 1 a) to indicate how much the central part of II Zw 40 is deviated from 
the reference radio-FIR relation defined for the BCDs. 



1.4 GHz flux from the entire galaxy (Fig. [3} indicates the presence 
of diffuse nonthermal emission in this galaxy. The trend that com- 
pact sources have lower r adio-to-FIR ratios is also found for lumi- 
nous infrared galaxies bv lCondon et alj dl99ll) . who also attribute 
the lower ratios to free-free absorption. 

To avoid free-free absorption, we use the 15 GHz luminosity. 
The best-fit model for t = 3 Myr is used for the 15 GHz luminosity 
of the II Zw 40 centre (Section 15. 2\ . Thi s lumi nosity fits the ob- 
servational data at 15 GHz bv lBecket al.1 ( 120021) . For com parison, 
we sho w the global luminosities of BCDs compiled by Hu ntet alj 
(2005) (Lfir is estimated by the method of lNagata et alj2002l ; see 
above). As we observe in Fig. |5p, there is a possibility that the 
radio-FIR relation is consistent with the global relation if the dust 
distribution is as compact as .Rdust = 100 pc (but note that the FIR 
luminosity is upper limit). If the dust distribution is as extended as 
-Rdust = 250 pc, the FIR luminosity is significantly suppressed (or 
the radio luminosity is enhanced) in the centre of II Zw 40. 

The deviation of the II Zw 40 centre in the ra dio-FIR rela- 
tion is actually expected, rather than peculiar. Indeed. lCannon et alj 
(2005) find a spatial variation of FIR-to-radio ratio by an or- 
der of magnitude within a metal poor dwarf galaxy, IC 2574. 
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iDumas et alj ( 120111) also show different radio-FIR relations be- 
tween spiral arms and interarm regions in M51. However, some 
results show that the radio-FIR correlation holds on 100 pc to 
sub-k pc scales JBeck & Goliall 19881: IxiTet al.lll992l : iHughes et all 
2006). There are some physical mechanisms that can contribute 
to the local variation of the radio -FIR relation such as dust pro- 
cessing by SN shocks ( ICannon et al. 2005), dust enrichment by the 
current star formation activity (Hirashita & Hunt 2008), amplifica- 
tion of magnetic fields feeck & Gollalll988l : I Dumas et al.luOl lh . 
diffusion and escape of cosmic ray elec trons jMurphv et alj|200q ; 
IHughes et al.ll2006r . lMurphv et alj(2008t) . etc. In particular, the dif- 
fusion of cosmic ray electro ns would predict a large FIR-to-radio 
ratio on small spatial scales JHughes et al. 2006). This is not likely 
to be the reason for the II Zw 40 centre though, as long as its stel- 
lar age is too young (< 3 Myr) for SNRs to generate cosmic ray 
electrons. 

If young active star-forming regions generally have similar 
radio-FIR emission properties to the central part of II Zw 40, free- 
free dominated submm emission and strong free-free absorption at 
centimetre and longer wavelengths will be a useful guide to find 
young active star formation in dense medium. This issue will be 
further investigated by collecting a sample of active star-forming 
BCDs with submm interferometry in the future. 



7 CONCLUSION 

In order to reveal the radiative properties of young active starburst, 
the central star-forming region in II Zw 40 was observed in the 
340 GHz (880 |J.m) band at ~ 5" resolution with SMA. A source 
associated with the central star-forming complex was detected with 
a flux of 13.6 ± 2.0 mjy and a 10-arcsec elongation in the east- 
west direction. The flux in the central part of II Zw 40 has been 
analyzed by usin g the theoretical radio SED model developed by 
iHirashita & Hung J2006I) . and interpreted along with interferomet- 
ric measurements at centimetre wavelengths in the literature. Then, 
we have found 

(i) that the SMA 880 (j.m flux is dominated by free-free emis- 
sion, and 

(ii) that possible contribution from dust emission to the SMA 
flux is less than 4 ± 2.5 mjy. 

Our models have been used to derive the radio-FIR relation 
of the II Zw 40 centre. We have discussed free-free absorption at 
low frequencies {v < several GHz; A > several cm) and spatial 
distribution of dust (i?dust) as possible factors affecting the radio- 
FIR relation. 
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